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ABSTRACT 

We present images of the young stellar object ASR 41 in the NGC 1333 star 
forming region at the wavelengths of Ha and [SII] and in the J, J, H, and if-bands. 
ASR 41 has the near-infrared morphology of an edge-on disk object, but appears 
an order of magnitude larger than typical systems of this kind. We also present 
detailed models of the scattering and radiative transfer in systems consisting of 
a young star surrounded by a proto-planetary disk, and the whole system being 
embedded in either an infalling envelope or a uniform molecular cloud. The best 
fit to the observed morphology can be achieved with a disk of ^200 AU diameter, 
immersed in a low density cloud {^2xl0~^^ gcm~^). The low cloud density is 
necessary to stay below the sub-mm flux upper limits and to preserve the shadow 
cast by the disk via single scattering. The results demonstrate that ASR 41 is 
probably not inherently different from typical edge-on disk objects, and that its 
large apparent size is due to the shadow of a much smaller disk being projected 
into the surrounding dusty molecular material. 

Subject headings: stars: pre-main- sequence — stars: formation — infrared: ISM 
— reflection nebulae 
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1. Introduction 

Young low-mass stars form by accretion of material through a disk, a process accompa- 
nied by outflow activity. When most of the nascent envelope has been dispersed in the later 
phases of star formation, a bipolar reflection nebula can be seen in the near infrared if the 
object is observed in the plane of the disk. Under this condition the star is obscured from 
direct view and therefore the total near-infrared flux is dominated by the scattered light 
from material above and below the disk plane. Many examples of such edge-on disks have 
been found recently. Most of the disks in those objects have radii inferred from the optical 
or near-infrared morphology that are of the order of magnitude of our solar system (i. e., 
^ 100 AU). The prototypical HH 30 IRS (Burrows et al. 1996) has a radius of ^ 250 AU, 
IRAS 04302+2247 (Lucas & Roche 1998) has a measured radius of ^ 200 AU in the near 
infrared, while the secondary component in the HK Tauri system (Stapelfeldt et al. 1998) has 
105 AU radius. The objects in the Taurus star forming region (Padgett et al. 1999) range 
in radius from 300 to 500 AU and the object found by Monin & Bouvier (2000) has about 
60 AU radius. In the Ophiuchus star forming region Brandner et al. (2000) found f«150AU 
disk radii, while Grosso et al. (2003) report a 300 AU radius and LkHQ;263C (Jayawardhana 
et al. 2002) has ^ 150 AU radius. 

Increasingly sophisticated models for these objects have been developed that treat mul- 
tiple scattering, absorption and radiative transfer by dust grains. These models reproduce 
many of the observed features (Whitney & Hartmann 1992, 1993; Wood et al. 1998, 2001; 
Whitney et al. 2003; Wolf, Padgett, & Stapelfeldt 2003), constrain the density distribution 
around the central star, and give information about the inclination of the disk. 

In the course of a near-infrared study of the star forming region NGC 1333 in the 
Per 0B2 complex, we found the young star ASR 41 (Aspin, Sandell, & Russel 1994) to show 
the morphology of an edge-on disk object with two exceptionally large scattering lobes. The 
edge-on disk structure of ASR 41 has been independently discovered by Elston et al. (2003) 
as part of their survey of giant molecular clouds. We assume NGC 1333 to be at the same 
distance of 316 pc that Herbig (1998) determined for the IC 348 cluster since they are both 
located in the same Per 0B2 molecular cloud complex. 

In this paper we present optical and infrared images of ASR 41 in combination with 
advanced scattering and radiative transfer models, and we argue that the unusual size of this 
object (radius ~ 3000 AU) is not indicative of the true extent of the disk, but rather is a 
projection effect of the shadow of a much smaller disk into the surrounding dusty molecular 
cloud material. 
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2. Observations and Data Reduction 

The NGC 1333 region was imaged with the QUIRC camera (Hodapp et ah 1996) at the 
UH 2.2m telescope, on the nights of Jan. 8 and 9, 2003 (UT), with the intent of surveying 
this region for substellar mass objects. The night was photometric, and observations in the 
J, H, and K "Mauna Kea Observatories" (MKO) filter set (Tokunaga, Simons, & Vacca 
2002) were obtained. The frames containing ASR 41 were part of a larger mosaic. The area 
surrounding ASR 41 was covered by 30 individual frames in each filter, the integration times 
were 150 s in J, 100 s in i/, and 50 s in K, for total integration times of 4500 s in J, 3000 s 
in H, and 1500 s in K. Photometric calibration is based on the UKIRT standard FSlll. A 
difficulty is posed by the bright sources and extended emission in this and all neighboring 
fields. The median filtering procedure to establish the sky frames leaves small residuals at 
the location of the bright sources, and contains large-scale gradients due to the presence 
of extended emission. The faint outer wings of the extended flux surrounding ASR 41 are 
therefore subject to artifacts from poor sky subtraction. 

The Ha (6569A/80A), [SlI] (6730A/80A), and /-band (8220A/1930A) images were ob- 
tained Oct. 13 and 14, 2001, using the 4 m Mayall telescope at NOAO with the MOSAIC 
CCD camera at a scale 0.26" /pix (Getman et al. 2002). /-band photometry of ASR 41 is 
based on the /-band magnitudes of neighboring stars given by Getman et al. (2002). Images 
taken earlier (Oct. 29, 1997) with the same equipment were used to measure the proper 
motion of the newly discovered Herbig-Haro object HH 727. 

3. Results and Discussion 

Fig. 1 shows the /i'-band image of ASR 41 and its surroundings, including the Herbig- 
Haro object HH 727, while Fig. 2 shows ASR 41 (and HH 727 as an insert) in the Ha, 
[SIX], /, J, //, and K bands. Going from the shortest wavelength covered here (Ha) to 
the J band, the morphology of ASR 41 changes dramatically. At Ha ASR 41 appears as 
a faint patch of nebulosity with a faint condensation at its center. At [SII], the central 
source appears brighter and the extended emission begins to show the "hourglass" shape of 
a bipolar nebula. The / band begins to show the bifurcation of the extended emission, albeit 
less pronounced than at the longer near-infrared wavelengths. The dark band at an angle of 
about 137°, bifurcating the reflection nebulosity associated with ASR 41, is most prominent 
from the J to the K band. The reflection nebula can be traced out to about 10" from the 
center (r^ 3000 AU) before it gets confused with artifacts from the sky subtraction process. 
The hnear size of the scattering region therefore is about 10 to 20 times larger than in the 
typical edge-on disk systems listed in the Introduction. Either, ASR 41 is a unique object 
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sur rounded by a huge extended disk or the dark band is the shadow of a much smaller (and 
typical) disk which is projected into the dusty material surrounding the object, as is strongly 
suggested by the wavelength dependence of its morphology. If the dark band were caused 
in its full extent by absorption in an edge-on dust disk in front of the scattering region, 
one would expect the dark bifurcating band to be more pronounced at shorter wavelengths. 
This is clearly not the case. Rather, our images suggest that the dark band is the shadow 
of a much smaller disk, projected into the surrounding dusty medium. The transition from 
single scattering at longer wavelengths, which preserves the shadow, to multiple scattering 
at shorter wavelengths, which fills in the shadow, occurs at wavelengths around 1/um. 

Within the dark band, a central object is clearly visible at all wavelengths from [Sll] 
to the K band. This central object has a FWHM of 0.69" in the K band, while the average 
FWHM of five other stars in Fig. 1 is 0.61±0.01". The central object is thus marginally 
different from a star; reflection nebulosity is contributing to its FWHM. Its position was 
measured in the J-band relative to 4 stars in the USNO-B catalog (Monet et al. 2003) to 
be a'* 28™ 5P.3 +31° 17' 40" (J2000), with estimated errors of f« 0.25" for the absolute 
coordinates. Astrometry of the central object in the [SII] image, using a larger set of USNO- 
B catalog stars than in the case of the infrared image, showed that it has the same coordinates 
as the corresponding object in the J-band image within the relative errors of ~ ±0.14". We 
don't see evidence for the wavelength dependent position of the central brightness peak 
that is often found in cometary or bipolar nebulae associated with Class I objects that are 
dominated by the absorption by a large and massive disk, e.g., (Hodapp et al. 1988). 

The central source has magnitudes in an 1.5" x 1.5" box centered on the central object 
of 7= 20.2, J= 17.4, iJ= 16.1, 15.3. The central source has a color of H-K = 0.73 and 
J-H — 1.36 while the color of the reflection nebulosity is bluer: We have taken measurements 
in 1.5" X 1.5" boxes at two additional positions to the north-east of the central star, i. e., 
in the reflection nebula, as indicated in Fig. 2. In the additional box closer to the central 
object, we measured H-K = 0.24 and J-H = 1.05, in the box farther to the north east we 
measure H-K = 0.34 and J-H = 1.15, with estimated errors in these colors of ^ ±0.05 mags. 
This is consistent with the flux from the central star and possibly the innermost parts of 
the reflection nebula being reddened by absorption in the small disk, and the flux in the 
extended reflection nebula being bluer, possibly because of preferential scattering by small 
dust particles. 

ASR 41 is not associated with a strong sub-mm source. The object lies just at the 
southern tip of a ridge of extended 850//m emission associated with the HH 12 complex 
(Sandell & Knee 2001), but it is not detectable as an individual point source in their map. 
From the fact that even the lowest contour (75 mJy/beam) of their map is not significantly 
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distorted at the position of ASR 41, we estimate an upper limit to its 850/xm flux of 200 
mJy. This imphes that the total mass of dust involved in the scattering of light from the 
central source, and the mass of the disk around the central star, must be relatively small. 
This point will be discussed in more detail in Section 4. 

To the south-west of ASR 41, a faint extended object at 2,^ 28"* 48M +31° 17' 10" (J2000) 
is most prominent in Ha and [SII] and faintly visible in K (Fig. 2. inserts). This is the pho- 
tometric signature of a low excitation Herbig-Haro object without significant [Fell] emission 
in the i/ band (Reipurth & Bally 2001). The object has therefore been named HH 727. 
Proper motion measurements on the [SII] images from 1997 and 2001 show a small (about 
1 pixel) shift of the photocenter of the HH knot towards the south west (PA = 215 ± 15). 
The shift corresponds to 100 km/s in the general direction away from ASR 41, roughly 
perpendicular to the plane of the bifurcating disk. While more precise measurements are 
clearly desirable, this strongly suggests that the HH object is physically associated with 
ASR 41. 



4. Modeling of the ASR 41 Disk Shadow 

We have computed detailed models for ASR 41 using a Monte Carlo scattered light code 
(Whitney & Hartmann 1992), (Whitney & Hartmann 1993); updated in (Whitney et al. 
2003). The models contain a number of simplifying assumptions, in particular about the 
grain properties, and therefore cannot be expected to reproduce all features of the observa- 
tions. The main purpose of the modeling is to demonstrate that the object can be explained 
well within the current paradigm of star formation. The models discussed here are based on 
the assumption that ASR 41 is essentially a Class II T-Tauri star surrounded by a disk of 
roughly the size of our solar system (r k, 100 AU), i.e., very similar to other edge-on disk 
systems. 

The two models differ in the assumptions about the distribution of the scattering dust 
surrounding the central star and its protoplanetary disk. Two types of models were consid- 
ered: disk plus constant density cloud (the "disk -|- cloud model") and disk plus infalling 
envelope (the "disk -|- envelope model"). Both reproduce the dark band as the shadow of a 
smaller disk. 

The geometry usually used to model young bipolar nebulae is an infalling envelope 
model, with the infall rate adjusted to model objects in different evolutionary classes (Whit- 
ney & Hartmann 1993; Whitney et al. 2003). For ASR 41, the infall rate is constraint by 
the upper hmit for the sub-mm flux. The "disk -|- envelope model" in Fig. 3 uses an outer 
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disk radius of 200 AU, a disk mass of 0.06 M©, an outer envelope radius of 7000 AU, and an 
infall rate of 3.5xlO~^M0/yr that is similar to the rate used in other successful models of 
young stellar objects (Whitney & Hartmann 1993). The infalling envelope has the density 
distribution of the free-fall models of Ulrich (1976) and Cassen & Moosman (1981). This 
model reproduces the basic features of the observations, but shows a steeper intensity gradi- 
ent with distance from the illuminating object than was observed. It predicts a 850/xm flux 
of 71 mJy, consistent with the sub-mm upper limit. 

The other model is based on the assumption outlined in qualitative form earlier that 
ASR 41 consists of a typical small cdgc-on disk object and a surrounding scattering cloud 
of uniform density into which the disk shadow gets projected. For this model, the density 
of the surrounding scattering material was adjusted so that single scattering dominates, 
which preserves the shadow effect and also coincidentally keeps the sub-mm flux within the 
observed upper limit. The "cloud model" uses a disk radius of 100 AU and a mass of 0.005 
Mq. for the best flt to the if-band data, the cloud itself is assumed to have a radius of 10000 
AU. Our "disk + cloud model" assumes a cloud density of 2x10"^*^ gcm~^, that is within the 
range of typical densities in large molecular clouds. This model produces a 850^m flux of 
40m Jy, well below the sub-mm flux limit. Increasing both the cloud density and disk mass 
by an order of magnitude leads to slightly better flts to the near-infrared fluxes, but predicts 
a 850/im flux of 234mJy, i. e., above the flux limit. 

Fig. 3 shows the iT-band image of ASR 41 (top, rotated to have the polar axis vertical for 
comparison with the models), and the two model flux distributions. Some noise was added 
to the model images to make them more visually comparable to the real image. The "disk 
-|- envelope" model is successful at matching the photometry, but morphologically, the flux 
is falling off too rapidly with distance from the central source due to the density distribution 
in the envelope. The "disk -|- cloud model" matches the near-infrared morphology better, 
including the large extent of the reflection nebulosity. 

Both model flts to the observed photometry indicate inclinations of about 80°, even 
though the match to the observed morphology, in particular the difference in flux from the two 
lobes could probably be improved by assuming somewhat lower inclinations. These problem 
our models have with precisely matching the near-infrared photometry and morphology of 
the object point to shortcomings in the dust model used here. Our code is currently limited 
to one dust model throughout, and given the extent of the object an interstellar grain model 
was chosen. This is a good approximation for the particles expected in a thin molecular 
cloud, but may not be a good model for the grains in a much denser, proto-planetary dust 
disk in the close vicinity of a young star, e.g. (Beckwith et al. 1990; Wood et al. 2002). Due 
to the obvious simpliflcations in our model, they are not suitable to derive the luminosity or 
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detailed information about the evolutionary state of the star at the center of ASR 41. 



5. Summciry 

We have presented images of the young stellar object ASR 41 in NGC 1333 at wave- 
lengths ranging from Ha to the if-band. The bipolar nebula seen at the longer of these 
wavelengths is much larger than typical edge-on disk systems but can be understood as 
the shadow of a smaller disk being projected into the dusty material of the surrounding 
molecular cloud. Detailed model calculations of the scattering and radiative transfer in this 
object have been presented. These models show that the morphological features and the 
near-infrared fluxes can be modeled by a small, low-mass disk and surrounding scattering 
medium without violating the upper limit on the sub-mm flux. ASR 41 is, most likely, a 
Class II T Tauri star with a disk of roughly the size of our solar system, and overall similar 
to other edge-on-disk objects. It is distinguished from those objects by the scattering of the 
disk shadow in the surrounding dusty molecular material. 

NOAO is operated by the Association of Universities for Research in Astronomy (AURA), 
Inc. under cooperative agreement with the National Science Foundation. 
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Fig. 1. — K-hand image of ASR 41. The objects are labeled by ASR number (Aspin, Sandell, 
& Russel 1994). The position of the Herbig-Haro object 727 is labeled HH. 
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Fig. 2. — Ha, [SII], /, J, H, and iT-band images of ASR 41. Inserted in each image is a 
smaller image of the HH-object 727, at the same spatial scale, but with the intensity scale 
stretched by a factor of two. The observed morphology of ASR 41 changes dramatically from 
Ha and [SII], where an extended nebula with a central, almost pointlike source is visible, 
to the / band, where the bifurcating dark band becomes noticeable. The bipolar nature of 
the object and the full extent of the dark band are fully established in the J band, and the 
morphology then remains virtually identical in the H and K bands. The inserted image of 
the HH object 727 is strongest in Ha and [SII], and indicated in K (due to H2 line emission), 
but is not detectable or very weak in the /, J, and H bands. Details of the morphology and 
photometry are discussed in Section 3. 
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Fig. 3. — Comparison of the K'-band image (rotated) of ASR 41 (top) with a "disk + envelope 
model" (center), and a "disk + cloud" model (bottom). 



